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We combine a light AsqTad antiquark with a nonrelativistic heavy quark to compute the decay constants 
of heavy-light pseudoscalar mesons using the ensemble of 3-flavor gauge field configurations generated by the 
MILC collaboration. Preliminary results for /s^ and fn, are given and status of the chiral extrapolation to /s 
is reported. We also touch upon results of the perturbative calculation which matches matrix elements in the 
effective theory to the full theory at 1-loop order. 



1. INTRODUCTION 

The "AsqTad" improved staggered quark ac- 
tion has made feasible the simulation of QCD 
with 3 flavors of sea quarks, with 2 of the fla- 
vors varying in mass from to^ to below ms/5. 
Unquenched simulation with this action at light 
sea quark masses produces agreement with exper- 
iment for several quantities: the T spectrum, B 
masses, and tt and K decay constants. Using a 
single set of input parameters, these quantities 
could not previously reproduce experiment [1]. 

Having removed this discrepancy and con- 
structed actions with good scaling properties, the 
uncertainties arising from chiral extrapolations 
can be studied more accurately. 

In recent work [2] we proposed and tested 
the use of improved staggered quarks as the 
light quark in heavy-light mesons. Since the 
heavy quark does not have the doubling problem, 
taste-breaking effects present in light staggered 
hadrons are suppressed in heavy-light mesons. 
Consequently the same operators used in Wilson 
fcrmion simulations can be used here, employing 
the identity between naive and staggered quark 
propagators. 

This work utilizes a subset of the public MILC 
configurations, the parameters of which we sum- 
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marize in Table 1; further details appear in [4]. 
The correct experimental kaon mass is obtained 
by tuning the valence strange quark mass to 
auQUiq = 0.040(1) on the au^nf^^ = 0.01 lattice. 
For clarity, we will quote quark masses in units 
of the mass, ms, corresponding to the physical 
strange mass. 

The NRQCD (improved through l/M^) and 
improved staggered actions are exactly as in [2] 
(see references therein). The operator matching 
has been carried out through 0{1/Mq) at 1-loop 
order in perturbation theory [5,6]. 

2. fs, AND 

The Bs decay constant is the simplest for us to 
compute since no extrapolations in valence quark 
mass are necessary and dynamical quark mass de- 
pendence is found to be small (see discussion be- 
low). The strange and bottom quark masses are 
fixed by tuning the bare masses to the physical K 
and T masses. On the lattice with rrf^l^ ~ TOs/4 
we find 

/s, = 260 ± 7 ± 26 ± 8 ± 5 MeV . (1) 

The first uncertainty is the statistical errors in the 
matrix elements and lattice spacing. The domi- 
nant uncertainty is the estimate of 0{al) effects 
neglected in the 1-loop matching calculation. The 
last two uncertainties estimate the errors due to 



Table 1 

Simulation parameters and inverse lattice spacing from two T splittings [3]. The lattice volume is 20^^ x 64. 

aupm^^^ quomf^ iVconf a-^2S - IS) g-^jlP-lS) 

0.01 0.05 568 1.59(2) GeV 1.58(3) GeV 

0.02 0.05 468 1.61(2) GeV 1.64(2) GeV 
0.03 005 564 1.60(3) GeV 1.68(4) GeV 
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Figure 1. Jh^ V'^ffs vs. inverse meson mass. 
Squares correspond to the light sea quark mass 
™u(? ~ "tns/i and diamonds to « TOs/2. The 
curve is the fit (to the squares) described in the 
text and the bursts are the fit results at physical 
values of mn ■ 



higher orders in the heavy quark expansion and 
discretization errors. 

Figure 1 shows the mass dependence of the 
combination fn, y/rriH^. We would like to extrap- 
olate to the charm region, where we are unable 
to calculate using NRQCD at this lattice spac- 
ing. Fits to a power series 
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give acceptable 's for > 2 - a linear fit gives 
an unacceptably large x^/DoF = 3. The extrap- 
olated value for fjj^ is 290 MeV. The largest un- 
certainty is again 0{a^) « 10% perturbative cor- 
rections. Other uncertainties are still being esti- 
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Figure 2. Light sea quark mass dependence on 
fn^ y/iTT-Hs ■ The heavy quark mass is « nib for 
the top 4 points and « 1.5 TOc for the bottom 4 
points. Octagons use from T(2S— IS) and 
diamonds use a"^ from T(1P — IS). 



mated. 

Within statistical errors, there is no sea quark 
mass dependence apparent in Fig. 1. In Fig. 2 we 
plot the same quantity for 2 values of heavy quark 
mass and with 2 definitions of the lattice spacing. 
Note from Table 1 that with m^^J" « (3/4)ms we 
see the reappearance of a scale ambiguity. With 
m^l « 'ms/2 the 2 lattice spacings differ by 1.5a. 
This effect apparently masks any sea quark mass 
dependence, as can be seen in Fig. 2, consequently 
we conclude that a chiral extrapolation in m^^^ 
will be an effect smaller than the other quoted er- 
rors. We emphasize that within statistical errors 
of 1.3% no scale ambiguity exists for m'^^ « ?7is/4 
which is where our result (1) is taken. 

3. CHIRAL BEHAVIOR 

The main benefit of using staggered fermions is 
being able to simulate closer to the chiral limit. In 
order to make maximal use out of the gauge field 
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Figure 3. ^$ = fB.^/mB'J lB,./nT^ plotted as 
a function of valence light quark mass in units of 



The crosses have 



and are uncor- 



related, and the squares have fixed nf^^ = ms/4: 
so are correlated. 



configurations, several values of valence quark 
mass niq are used. So far we have accumulated 
data with nf^l^ > nis/A and nuq > m^/S, com- 
pared to TTiq > 0.7 TTis which is the state-of-the-art 
using Wilson- like fermions [7] . 

In Fig. 3 we plot fsjlsg (times ^/m^Jm^ 
or 1.01) against vriq/ms. Crosses are unquenched, 
except that the dynamical strange quark mass 
is slightly heavier [rrf^^ = (5/4)tos) than the 
valence strange quark mass (m™' = rus), and 
squares have light sea quark mass fixed to ms/4. 
A partially quenched analysis will utilize the cor- 
relations between data points computed on the 
same configurations. 

Becirevic et al. recently noticed that the coef- 
ficient of the chiral logs in the double ratio 
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is numerically smaller than in either ratio individ- 
ually [8]. Combining our results with the MILC 
collaboration results [9] for /at/Ztt (which is 1.22 
experimentally) yields Fig. 4. Extrapolating R to 
r = 0, a range of 1.0 - 1.1 would correspond to 
^ between 1.22 and 1.34. [^BbJBb^ = 1.01(3) 
[7].) Much work remains to be done before we 
have a trustworthy and precise final result, but 
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Figure 4. The double ratio R (see Eq. (3)) plotted 
as a function of valence light quark mass in units 
of rUs- Symbols are as in the previous figure. 



these simulations with masses below ms/2 should 
cast new light on the chiral extrapolation of this 
important quantity. 
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